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Abstract 

Recent experimental results suggest that metallic liquids universally exhibit a high-temperature 
dynamical crossover, which is correlated with the glass transition temperature (Tg). We demon¬ 
strate, using molecular dynamics results for Cu64Zr36, that this temperature, Ta ^ 2 xTg, is linked 
with cooperative atomic rearrangements that produce domains of connected icosahedra. Super¬ 
cooling to a new characteristic temperature, To, is shown to produce higher-order cooperative 
rearrangements amongst connected icosahedra, which manifests as the formation of large Zr-rich 
connected domains that possess macroscopic proportions of the liquid’s icosahedra. This coincides 
with the decoupling of atomic diffusivities, large-scale domain fluctuations, and the onset of glassy 
dynamics in the liquid. These extensive domains then abruptly stabilize above Tg and eventually 
percolate before the glass is formed. All characteristic temperatures {Ta, Td and Tg) are thus 
connected by successive manifestations of the structural cooperativity that begins at Ta. 

PACS numbers: 
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I. INTRODUCTION 


As a liquid is supercooled below its melting temperature, the characteristic timescales 
of its dynamics become increasingly stretched. In a dramatic departure from exhibiting 
simple diffusive particle dynamics, its viscosity, rj, increases rapidly as it develops general, 
glassy dynamic features^^-. Efforts to universally describe these phenomena as a liquid is 
supercooled to its glass transition temperature {Tg) remain divided in both their approaches 
and successes^"—. In these theories, the roles played by structure, if any, are a major 
point of discrepancy. The molecular dynamics (MD) results presented here reveal a striking 
correlation between structure and dynamics across a broad temperature range. In particular, 
a high temperature structural crossover is shown to underlay a recently reported dynamical 
one^^. 

Numerous works point to the fundamental importance of locally-preferred structures 
(LPS)s in glass forming liquidsl^'i^^— . One such thermodynamic theory describes the growth 
of frustration-limited domains (FLD)s of the EPS as the liquid is cooled^'i^.'^. This FED 
theory accurately describes glassy relaxation processes and predicts a temperature, T^, 
below which a fragile liquid develops a super-Arrhenius temperature dependence for rj. Very 
recently, an analysis of experimental data for a wide variety of metallic liquids revealed that 
Ta correlates strongly with Tg, such that ^ 2 x Tg—. This temperature may correspond 
to a crossover phenomenon created by competition between conhgurational excitations and 
phonons^i. Remarkably, it was found that, by using Ta and particle density as scaling 
factors, the viscosity data for all strong and fragile metallic liquids fall on a universal curve 
across a broad temperature range^^l^. 

Despite the merits of the FED theories, an important problem persists. The measured 
structural changes of the liquid across the accessible supercooled temperature range are ex¬ 
tremely subtle, casting doubts on the extent and even the existence of the FEDs2i>^ as well as 
the structural changes that accompany processes at Ta^'^. Therefore, observing FEDs and 
studying their role in the dynamical crossover at Ta would bolster FED theory, the inter¬ 
pretation of the apparent universality across metallic liquids, and the overall understanding 
of the relationship between structure and dynamics in liquids 

In this work, we demonstrate the development and growth of FEDs in a metallic liquid 
as it is cooled through Ta, providing a coherent picture of the structural and dynamical 
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features of the cooperativity that arises at this crossover temperature as well as within 
the supercooled regime. Rapidly growing fluctuations in domain sizes are linked to the 
development of glassy dynamics, leading to the identihcation of another crossover at lower 
temperature, Tp, which may be possible to observe in experimental structural studiesSSi^. 
While supercooling, these extensive domain then quickly stabilize, marked by a sudden drop 
in domain fluctuations, and percolate across the system before the temperature reaches Tg. 
The MD results presented here suggest that all of these characteristic temperatures, Ta, To 
and Tg correspond to a cascade of cooperative structural rearrangements involving a LPS 
that begins at T^, a result consistent with the prediction from recent experimental studies 
of viscosity^. 

The remainder of this paper is organized as follows: in section II, we introduce the meth¬ 
ods of our MD simulations. In section III, we discuss our methods of our Voronoi analysis, 
and introduce dehnitions of icosahedral ordering. Section IV introduces the temperatures 
Ta and Td, and the dynamical features that accompany them. Section V dehnes liquid re¬ 
laxation and restructuring timescales and discusses the role played by cooperative structural 
rearrangements at Ta and Tp. Here also we show that Ta marks the onset of connected 
icosahedral ordering. Section VI links that dynamical features that begin at with the 
rapid proliferation of connected domains of icosahedra. Section VII provides an overview 
of the structure-dynamic relationships of the liquid and its icosahedron domains from high 
temperatures down to Tg, where were hnd that a single domain of icosahedra has percolated 
the system. Section VIII contains a summary of our hndings. 

II. METHODS 

Results were obtained from classical MD simulations of the metallic liquid Cu64Zr36, 
which among best glass forming compositions of this alloy^. All MD simulations were 
conducted using LAMMPS^ with GPU-parallelization packages^. Atomic interactions in 
the Cu64Zr36 system, including approximations of many-body interactions, are described 
using the Finnis-Sinclair generalization of the embedded atom method^I. This semi empirical 
potential was selected for its ability to accurately reproduce structures and properties of 
Cu64Zr36 in both the liquid and glass regimes, as conhrmed by comparisons made with x-ray 
diffraction data and ab initio calculation results^. The reported calorimetric glass formation 
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temperature {Tg) for this potential is approximately TSOii'. We confirm this by measuring 
the system volume, given zero average pressure, as a function of temperature, and identifying 
the temperature at which the thermal expansion coefficient changes. 

Isothermal, isobaric simulations (sampled from ensembles with constant iV, (P), and 
(T)^) were conducted with (P) = 0 and N = 3 x atoms, using periodic boundary 
conditions, the Verlet integration method, and the Nose-Hoover thermostatting method. 
Each simulation was initialized with a random atomic configuration that evolves at T = 
3300P for 0.5ns to achieve equilibrium. The liquid was then quenched at a rate of lO^^P/s 
down to its target temperature and subsequently relaxed for 20ns before structure and 
dynamics data were collected. The integration timestep during the quenching and relaxation 
processes was 5/s. The integration step was then decreased to 2/s, and data was collected 
for each timestep. For a given target temperature, atomic-level data, such as the positions 
and velocities were recorded across a 0.2ns interval; the system’s Cauchy stress tensor was 
recorded for 8ns. Additional independent simulations were conducted for all temperatures 
below llOOiC in addition to 2200P and 2500P. To check for system size effects, additional 
simulations of = 1 x 10^ atoms were performed. 

III. ICOSAHEDRAL ORDERING AND FRUSTRATION-LIMITED DOMAINS 

Voronoi analysis utilizing a weighted-bisector method was performed utilizing the Voro-|--|- 
software library^. Small faces comprising less than 0.5% of a Voronoi cell’s surface area were 
removed. The analysis was performed on each time step for which atomic positions were 
recorded. The LPS for Cu 64 Zr 36 is an icosahedron-shaped cluster of atoms^“— , consisting 
of a Cu atom with 12 nearest neighbor atoms, which form a weighted-Voronoi cell with 
12 pentagonal faces. Distorted icosahedra, which contain some non-pentagonal Voronoi 
faces, are treated as distinct structures and are not considered to be LPSs. FLDs are 
thus comprised of connected icosahedra. Two icosahedra are connected if they share a 
vertex, an edge, a face, or they interpenetrate, meaning that they share 1, 2, 3, or 7 atoms, 
respectively^^. 
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FIG. 1: (Color online) Viscosity, rj, as a function of y. r/ begins to exhibit super-Arrhenius growth 
below Ta = IbSOA'. (Inset) The percent residual of rj from the Arrhenius fit. 

IV. DYNAMICAL FEATURES OF LIQUID CuuZrs^i: BREAK DOWN OF 
STOKES-EINSTEIN AND DECOUPLING OF DIFFUSIVITIES 

Before analyzing the strnctnral evolntion of Cn 64 Zr 36 , a snmmary of the dynamical fea- 
tnres that develop in the liqnid across a broad temperatnre range is presented. Viscosity 
[r]) data obtained nsing the Green-Knbo relationship for the stress tensor antocorrelation 
fnnctioii^ are shown in the main panel of Fignre[Il A change from an Arrhenins tempera¬ 
tnre dependence to a snper-Arrhenius form occnrs when the liqnid is cooled below 1550K, 
which is therefore identihed as Ta- In accord with recent experimental study^, this valne is 
near 2 x Tg {Tg ^ TSOK— ). 

The liqnid diffnsivity data, plotted in the inset of Fignre|2], maintains an Arrhenins form 
throngh lower temperatnres than does r], indicating that Ta also marks the breakdown of 
the Stokes-Einstein relationship for Zr atoms, 

V-DzrC^S -, ( 1 ) 

Kzr 

which are the larger diffnsing solnte species in the system^>^. Here, Dzr is the diffnsion 
coefficient of tagged Zr particles as measnred from the long-time integral of their velocity 
antocorrelation fnnctions, T is the temperatnre, and Rzr is the effective particle radins. The 
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FIG. 2: (Color online) The Stokes-Einstein ratio as a function of T. The Stokes-Einstein rela¬ 
tionship becomes violated below (Inset) Cu and Zr diffusivities versns A transition from a 
high-T linear fit on a log scale (solid line) occnrs at To = lObOitT for both species. 

data presented in the main panel of Fignre [2]shows that for high temperatnres the Stokes- 
Einstein ratio is ronghly a constant 0.67, which is comparable to ~ 0.65 1 /At T/i, 
the accelerated growth in rj causes a deviation from this constant, indicating that the liquid 
is departing from its “simple” behavior, and is developing dynamical heterogeneities^'^. 
Specihcally, slow-moving regions develop in the liquid that dictate the timescale of structural 
relaxations, as measured by rj, whereas fast-moving regions control the dynamical timescale 
proportional to A. 

It is striking to hnd that the Stokes-Einstein relationship is violated at such a high 
temperature - literature reviews of supercooled liquids broadly indicate that this relationship 
is expected to hold well into a liquid’s supercooled regime^*^. However, this prescription is 
chiefly informed by empirical studies of molecular liquids, such as o-terphenyl, where the 
violation occurs near 1.2 x Tg (near the mode-couple temperature)^!^. Our Ending, that 
the Stokes-Einstein relationship is violated well above 1.2 x Tg, are supported by several 
empirical studies of metallic liquids^^— . Furthermore, an illuminating concurrent MD study 
of liquid Cu 4 oZr 5 iAl 9 submitted very shortly after the original submission of the current work 
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FIG. 3: (Color online) (Top) The mean-squared displacement time-trajectories of atoms for various 
liquid temperatures on a log-log scale. The temperatures included are separated by lOOii' intervals. 
A distinct plateau feature, characteristic of glassy dynamics, develops for T < Tq. (Bottom) The 
ratio of the diffusivity of Zr atoms to that of Cu atoms. The diffusivities are found to decouple 
sharply for T < Tq. 

similarly found that the relation breaks down specihcally at ^ 2 x Tg for that alloy as 
well^. Taken together, these results bolster our finding of Stokes-Einstein violation at Ta 
to a broader range of materials 

It is important to note that the local violation of the Stokes-Einstein equation at Ta 
does not signify the onset of strong glassy dynamics in Cu 64 Zr 36 . That is, the system does 
not yet exhibit plateau-separated fast {/3) and slow (a) relaxation processes at Ta- This 
can be seen in the inset of Figure [2l which displays the temperature dependence of the 
diffusivities for Cu and Zr particles, and in the top panel of Figure [HI which contains a 
series of mean-squared displacement (MSD) time trajectories for several temperatures. The 
MSD trajectory at Ta exhibits an immediate transition between simple ballistic and diffusive 
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motion regimes without any plateau feature. This changes once the diffusivity, D(T), departs 
from an exponentially-decaying function of 1/T at lOSOi^, which we label Tq. D reflects the 
long-time asymptotic behaviors of the particles’ MSDs: (5r(t)^) ~ 6Dt {t —)■ cxo), thus the 
accelerated decline in D(T) for both atomic species reflects the development of a sustained 
plateau that separates two-step relaxation processes in the supercooled liquid (Fig[2]), as well 
as growing diffusive-motion activation energy barriers. This is characteristic of caged particle 
dynamics and is an important equilibrium signature of the impending glass transition!^>^. 

Also at Td, the motions of Cu and Zr atoms begin to decouple, leading to a sharp decrease 
in with decreasing temperature; this is shown in the lower panel of Figure [31 Here, the 
system’s energy landscape begins to to disproportionately impede the activation of the larger 
Zr atoms, whereas above To, the two species diffuse via the same physical mechanism with 
^ -Rcvjs^ rppjg decoupling strongly suggests that a structural feature begins to emerge 

^Cu -t^Zr 

at To that skews the supercooled liquid’s dynamics. We will investigate this underlying 
structure in section VI. Lastly, we note that a recent empirical study of a Zr-based metallic 
liquid reports a similar decoupling of dynamics at a temperature comparable to To^- 
In summary, Ta (1550iF) is the temperature below which r] begins to exhibit super- 
Arrhenius growth, and where the Stokes-Einstein relationship hrst breaks down in the liq¬ 
uid. To (lOSOiF) then marks the early onset of stretched, glassy dynamical features in the 
supercooled liquid, which appear to be shaped by an emerging robust structural feature. 

V. Ta and cooperative structural rearrangements 

The crossover from the “simple liquid” regime at Ta marked by the development of super- 
Arrhenius relaxation times, and the development of glassy dynamics at To both coincide 
with the onset of structural cooperativity associated with FLDs of connected icosahedra. 
Our analysis of this cooperativity follows and expands on the work of Iwashita, Nicholson, 
and Egami^ and introduces three pertinent timescales. First, the Maxwell time, or shear 
stress relaxation time, 

provides a timescale during which the liquid behaves like a solid and exhibits an elastic 
shear response without flowing^>21. aij{t) is the time-dependent Cauchy stress tensor—. 
Second, the local cluster time, tlc, is dehned as the average time required for an atom 
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FIG. 4: (Color online) Ratio of the Maxwell time to the local cluster time and the domain-connected 
icosahedron lifetime, respectively, tm surpasses tlc and at Ta and To, respectively. (Inset) 

A 2D illustration of the cooperative rearrangements that begin contributing to the liquid relaxation 
process at Ta and at To, respectively. Here the 2D pentagons represent 3D icosahedral clusters of 
atoms. 

to change its coordination number, i.e. the time needed for an atom to lose or gain a 
Voronoi neighbor. In addition to these two times, we introduce the icosahedron lifetime, 
Ticos, which is the average time required for an icosahedron to lose a vertex atom, gain an 
extraneous vertex atom, or distort its shape. Tjcos for an icosahedron involved in a domain 
via interpenetrating connection is labelled ■ Measurements of tlc and Ticos provide 

a context for understanding the timescales over which atomic rearrangements of the liquid 
structure occur. 

The primary portion of Figure 0] shows the ratios of tm to tlc and respectively, 

as a function of temperature. Liquid relaxation in the high temperature, simple liquid 
regime (T > Ta) is characterized by uncorrelated local rearrangements of atomic clusters: 
Tm ~ Tlc- That the high temperature liquid relaxation time is generally controlled by tlc 
in atomic liquids was hrst shown by Iwashita et Here, tm is so short that consecutive 
atomic rearrangements occur only after the surrounding environment has relaxed in response 
to shear strains. It is not until the system is cooled through Ta that tm > tlc, and that 
typical atomic clusters can begin to rearrange with one another on a time scale during which 
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FIG. 5: (Color online) The populations of connected and isolated icosahedra, plus their sum. A 
low-T linear fit on the log scale is indicated with a dashed line. A crossover in the isolated and 
connected populations occurs as the liquid is cooled throug Ta- (Inset) The average connected- 
icosahedron domain size, s, on a log scale. A domain’s size is given by the number of icosahedra 
participating in a singly-connected structure 

the liquid behaves like a solid. Here cooperative restructuring can begin contributing to the 
liquid’s relaxation process^i. 

What is produced by these cooperative rearrangments? That is, what difference is there 
between the average structure of the liquid at temperatures above Ta and those below? 
Measuring the the populations of the isolated, and connected icosahedra in a volume con¬ 
taining 3 X 10^ atoms reveals that a structural crossover occurs at Ta'- above Ta it is more 
likely to find isolated icosahedra than connected ones, and below Ta connected icosahedra 
become the dominant population. This can be seen in Figure |5], which shows, on a log scale, 
the evolution of these icosahedron populations as well as the total number of icosahedra. 
The total population of icosahedra proceeds to grow exponentially with inverse temperature 
for Ta > T. In the context of our timescales, we find that in the regime where Ta > T > T^ 
icosahedra are stable in comparison to the average local cluster configuration (t^^^ > t^c), 
and are relatively robust on the timescale of liquid relaxation ^ tm) across this tem¬ 

perature range. Thus the liquid relaxation process chiefly involve non-icosahedral atomic 
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FIG. 6: (Color online) The fraction of the total number of icosahedra in a 30,000 atom system, 
^icosj contained in an FLD. Below Tq, a domain contains a macroscopic number of the system’s 
icosahedra. Tp^^k marks the maximum susceptibility of the fractional occupation. 

clusters, which cooperatively rearrange amidst the relatively inert icosahedra and form off 
of them, new, connected icosahedra. This process is depicted in a simplihed 2D illnstration 
in the inset of Fignre HI 

It is qnite remarkable that the distinct dynamical featnres associated with are linked 
with the incipient icosahedral ordering in this system. While the strnctnre-dynamics featnres 
of icosahedron networks are becoming well-docnmented for temperatures near Tg— the 
roles played by these structures at higher temperatures have remained largely unknown. 
Thus our hndings stand amongst the hrst explicit accounts of the structural features that 
accompany the emergence of solidlike features in high temperature metallics liquids. They 
also correspond closely with the narrative of FLD theory, which predicts the development of 
frustrated domains at Ta in association with an avoided critical point in the system^'I^'^. 
We proceed by further tracking this icosahedral ordering into the supercooled region, where 
prominent structural features are expected to emerge near To. 
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FIG. 7: The concentration of Cu atoms in: an average 13-atom icosahedron (red dashed-line), the 
connected icosahedron domains (red circles), and in non-icosahedron structures (yellow squares). 
When approaching Tg, the preferential sharing of Zr atoms amongst icosahedra allows the extensive 
FLDs to approach the overall alloy composition, while preserving the Zr-rich environment required 
by the icosahedra. 

VI. Td and the rapid onset of dominant icosahedral ordering 

As the liquid is supercooled through To, its relaxation time grows such that tm 3> 
for T <Td (see FigHj), and thus its relaxation processes begin to involve consecutive struc¬ 
tural rearrangements of connected icosahedra within FLDs. The illustration in Figure 0] 
depicts a simple schematic of such processes. We hnd that the onset of this medium-ranged 
cooperative restructuring accompanies the enhancement of the diffusive-motion energy bar¬ 
rier and the early development of general glassy dynamics observed at T^, and it leads to a 
dramatic proliferation of highly connected FLDs of icosahedra throughout the liquid. 

The size of an FLD is measured by the number of icosahedra participating in the con¬ 
nected domain. The inset of Figure [5] shows the the weighted average FLD size, s, whose 
growth accelerates tremendously with decreasing temperature - note that the data is plot¬ 
ted on a log scale. The fractional occupation of an average FLD - the proportion of the 
icosahedron population in a 30,000 atom system contained in a typical connected domain - 
is plotted in Figure [6] versus inverse temperature. Here we see that cooling through To leads 
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to a dramatic unification of the liqud’s FLDs. Above To a typical FLD contains a neglible 
portion of all the icosahedra, whereas below To the rapidly growing FLDs begin join with 
one another, such that macroscopic number of the system’s icosahedra can be found in a 
single connected FLD. The peak value in the susceptibility of the fractional occupation num¬ 
ber, -^{s/Nicos), resides near 900iF and is marked Tp^ak- By 825K, s/Nicos nearly reaches 
unity - only a negligible number of icosahedra can be found to be separate from a single, 
extensive network of icosahedra. 

These extensive FLDs, we propose, are the underlying structures responsible for the de¬ 
coupled dynamics and the enhanced barriers to diffusion that emerge below T^. Indeed, the 
disproportionately depressed diffusivity for Zr can be accounted for by considering the chem¬ 
ical ordering that accompanies the formation of FLDs. In accordance with efficient packing 
schemes of different-sized hard sphere8^>^, icosahedra are Zr-rich relative to the overall 
alloy composition Cu 64 Zr 36 - below an icosahedron has an approximate composition of 
Cu 57 Zr 43 . An earlier MD study showed that non-interpenetrating connections between the 
icosahedra in an FLD strongly favor sharing Zr atoms amongst one another—; for instance, 
80% of all shared vertices are Zr atoms, whereas the icosahedron shell composition is only 
43% Zr. Thus non-interpenetrating connections within FLDs are utilized to foster a locally 
Zr-rich environment for the icosahedra, while the FLDs at large can, necessarily, maintain 
a composition closer to that of the alloy. The temperature dependent Cu concentrations 
of individual icosahedra (including the center Cu atoms), FLDs, and the non-icosahedron 
regions of the liquid are plotted in Figure [71 

There are numerous studies that demonstrate the unique role played by icosahedron 
FLDs in creating strong dynamical heterogeneities in Cu 64 Zr 36 as it is supercooled towards 
Tg—^^. In particular, Y. Zhang et al. recently showed quite explicitly that diffusing atoms 
strongly avoid regions of icosahedral ordering, and that they instead reside in liquidlike re¬ 
gions of the system^. They demonstrate that underlying the development of medium-ranged 
icosahedral ordering and glassy dynamics is a crucial correspondence between dynamical and 
structural correlation length scales in the liquid. We thus conclude that the growing barrier 
to diffusion and decoupled motions of Cu and Zr atoms, which ensue below To, is created 
by the rapidly proliferating, Zr-rich FLDs of connected icosahedra that form amidst Cu-rich 
liquidlike regions in the system. 
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FIG. 8: Relative fluctuations of the mean connected-icosahedron domain size over time, as a 
function of scaled by y/N, where N is the number of atoms in the system. A deviation from 
the high temperature Arrhenius fit (dashed curve) begins at Ta- Domain-size fluctuations grow 
rapidly upon cooling through To, and become system size independent. The distinct peak near 
925Ar demonstrates the abrupt stabilization and growth of a large FLD, which percolates above 
Tg. Error bars are derived from averaging two independent runs 

VII. FLD FLUCTUATIONS AND FALLING OUT OF EQUILIBRIUM 

Both cooperative characteristic temperatures, Ta and Td, manifest as enhancements in 
fluctuations of the average domain size, s . Figure El shows the relative time-fluctuations 
in domain size, versus inverse temperature for system sizes of 1 x 10"^ and 3 x 10^ 
atoms. The fluctuations are scaled by a factor of \/N /lOO to make-clear the behavior of 
system size effects. Here, s{t) is the mean domain size of the system at a time t, and a-g 
is the standard deviation of s{t) from its time-averaged value, s. As the liquid is cooled 
through Ta, fluctuations in the domain size become enhanced, growing faster than does the 
high-temperature Arrhenius £t (dashed curve) with decreasing temperature. Above To, the 
scaled fluctuations for the two system sizes collapse onto a single curve, indicating that FLD 
fluctuations scale as ~ l/y/N in this temperature range. Cooling further through Td leads 
to a breakdown in this scaling law, with the peak fluctuations growing as ~ \/N on the 
scaled axis - that is, the peak-value fluctuations appear to become system size independent. 
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By T = 925K the cooperatively-rearranging FLDs attempt to grow far beyond their average 
sizes before inevitably collapsing, reaching relative fluctuations of 70%. 

This behavior changes abruptly as the liquid is supercooled through 92577. Here the 
FLDs begin to stabilize and the fluctuations in FLD size drop rapidly with decreasing 
temperature. We must consider whether this drop in fluctuations is actually a mark of 
structural stabilization or if it is a mere consequence of our simulation timescale becoming 
insufficient for capturing metasable equilibrium dynamics. It is certainly true that the 
fluctuations must inevitably fall to zero as we approach Tg given that tm will far outgrow 
our simulation timescale. That being said, tm is still sufficiently short at 90077, where the 
fluctuations have already started to fall, that our liquid has been allowed to evolve for a 
time four orders of magnitude longer than the liquid relaxation time, and our data collection 
window is two orders of magnitude longer. Furthermore, several studies show that utilizing 
slower quench rates, longer relaxation time, and thermalization techniques only further 
enhances the degree of icosahedral ordering and the connectivity of the resulting FLDs^^i^. 
This indicates that the drop in fluctuations indeed reflects pronounced stabilization of the 
FLDs, and that the fluctuations in the extensive FLD’s size would approach zero even if we 
were able to remain in equilibrium near Tg. By 90077, nearly 40% of the system’s icosahedra; 
as we saw in Figure [6l this fraction continues to grow rapidly until a single FLD percolates 
the system by 82577—, where the fluctuations nearly reach zero. A structure is said to 
percolate a material when a singly-connected domain of that structure spans the material 
(e.g. a connected domain of icosahedra is as long as the diagonal length of our simulation 
box). Below Tg, the percolated FLD shapes the properties, mechanical and others, of the 
glass. This has been demonstrated by a number of excellent studies for several glasses^i^'— . 
Not discussed in this paper are the roles played by the various types of connections that can 
form between icosahedra. For detailed discussions of this important aspect of icosahedral 
ordering, we refer the reader to early studies of this subject^'^. 

In the limit of large system size, the peaked fluctuations in Figure [H] resembles closely the 
susceptibility of the fractional occupation plotted in Figure El We note that this and other 
features of the rapid development of extensive icosahedral ordering in this system indicates 
the possible presence of a liquid-liquid phase transition above Tg—. This will be investigated 
in detail in a future publication. 

Lastly, the time-dependent behavior of the FLDs that produce the fluctuations reported 
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FIG. 9: (Color Online) The connected-icosahedron domain sizes present at a given simulation time 
step plotted versus time. System temperatures are indicated on the right side of the plots. On the 
Ips timescale, an orange trajectory traces the evolution of the largest domain size. 

in Figure [8] are shown in Figure |9l Here, we plot all of the domain sizes that are present at 
each simulation time step, for temperatures Ta, Td, Tp^ak, and Tg. The left side and right 
sides of the plot show this time dependence on Ips and lOOps timescales, respectively. In 
the left panel, the evolution of the largest FLD size is highlighted by an orange trajectory. 
With decreasing temperature, down to Tpeak, a clear enhancement in the time variation of 
FLD sizes is observed. At Tg, & single, large FLD dominates the glass and exhibits minimal 
fluctuations across the 200ps time window. 


VIII. SUMMARY 

In summary, we have employed MD simulations to show that distinct structural features 
accompany dynamical crossovers in liquid Cu 64 Zr 36 . Comparing the liquid relaxation time, 
tm, to measures of structure lifetimes, tlc and allows us to identify the temperature 

ranges in which cooperative structural rearrangements contribute to the liquid’s relaxation 
processes. At Ta ^ 2 x Tg atoms begin to cooperatively form connected icosahedra off of 


16 





























isolated ones, serving as the beginning of the icosahedral ordering that eventually dominates 
the system. This accompanies the onset of the super Arrhenius growth of rj and the break 
down of the Stokes-Einstein relationship in the liquid. 

In the supercooled regime, Tj:, ^ 1.4 x Tg marks the early development of stretched glassy 
dynamics, and the decoupling of Cu and Zr diffusivities. These dynamical features are a 
manifestation of rapidly growing Zr-rich domains of connected icosahedra, which possess 
macroscopics numbers of the liquid’s icosahedra. Below 925K these FLDs stabilize abruptly, 
and eventually percolate the system before the liquid reaches Tg. Our results provide needed 
evidence and explanation for the structural roles played at T^, which according to recent 
experimental results is an important characteristic temperature for all metallic liquids^^, 
and more broadly illustrates that a cascade of cooperative structural and dynamical effects 
begins at this temperature and characterizes the supercooled liquid as it approaches Tg. 
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